Phonon assisted dynamical Coulomb blockade in a thin suspended graphite sheet. 
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The differential conductance in a suspended few layered graphene sample is found to exhibit a 
series of quasi-periodic sharp dips as a function of bias at low temperature. We show that they can 
be understood within a simple model of dynamical Coulomb blockade where energy exchanges take 
place between the charge carriers transmitted trough the sample and a dissipative electromagnetic 
environment with a resonant phonon mode strongly coupled to the electrons. 
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One of the great challenges of molecular electronics is 
to access electron-phonon coupling at the single molecule 
level. Mechanically tunable atomic break junctions with 
trapped small molecules such as {H2, D2, H2O) have been 
shown to exhibit a spectroscopic signature of their char- 
acteristic phonon modes [ij. The signature of phonons 
is also spectacular in the Coulomb blockade regime for 
a molecular single electron transistor: The typical reso- 
nant tunneling peaks as a function of gate or source-drain 
voltage are surrounded by satellites, which correspond 
to the emission or absorption of one or several phonons. 
Specific vibrational modes were identified in this way in 
fuUerenes and suspended carbon nanotubes [2,l3|,y]. The- 
oretical models [5|, |6|, [Zl were developed to describe these 
vibrational side bands in molecular transistors, involving 
either a quantum or classical treatment of the electron 
phonon coupling. In all these investigations the single 
level spacing within the molecule is larger than the en- 
ergy of the vibration modes coupled to the molecule, so 
that only a single molecular level needs to be considered. 
In the present work we investigate the opposite limit of 
a mesoscopic dot where both the single level spacing and 
Coulomb charging energy are smaller than the energy 
of the phonon mode considered. Moreover the trans- 
mission of the electrodes corresponds to an intermedi- 
ate tunneling regime described by the physics of dynam- 
ical Coulomb blockade (DCB). The samples are micron 
size few atomic layer graphite foils suspended between 
two platinum electrodes. The differential conductance 
exhibits at low temperature around zero bias a power 
law increase characteristic of DCB through the contacts. 
The graphite foil sample itself constitutes the electro- 
magnetic environment. More original, on the thinest 
sample (with 30 graphene layers) a series of periodic 
replica of this Coulomb blockade anomaly was detected 
at multiples of 20 meV, corresponding approximately to 
the lowest energy out of plane optical mode in graphite 
(Z0')'8]. These sharp dips were not observed on two con- 
trol graphite samples which were likewise suspended, had 
similar resistance and lateral dimensions but were more 
than 30 times thicker. We analyse these results with a 
simple model, inspired by [7| , of a mesoscopic island con- 
nected to electrodes via tunnel barriers. We model the 



island by a continuous electronic spectrum coupled to a 
single phonon mode, which leads to an oscillating trans- 
mission of the barriers at the contacts [7|. This model 
can also be solved using the so called P{E) theory de- 
veloped by Ingold and Nazarov [9| to describe energy 
exchanges of a Coulomb blocked tunnel junction with a 
dissipative electromagnetic environment and presents a 
striking analogy with the behavior of a tunnel junction 
coupled to an electromagnetic resonator in series with an 
ohmic environment as investigated by [l0|. We finally 
deduce from the field dependence of these dips an order 
of magnitude for the relevant electron-phonon coupling 
parameter in the system. 



The samples were prepared by exfoliation of a highly 
oriented pyrolytic graphite (HOPG) single crystal and 
deposition across a slit etched in a silicon nitride mem- 
brane separating two Pt metallic contacts. The num- 
ber of graphene layers was estimated from transmission 
electron microscopy pictures, see figH] for the thinnest 
sample, which contains between 25 and 35 layers. The 
electrical contacts were obtained just by pressing the 
sample onto the electrodes. Thus the sample resistance, 
Rt — lOOfcil at 4.2 K, mainly consists of the resistance at 
the contact, and increases as the temperature is reduced. 
The square resistance of the graphite layer itself is not 
expected to be larger than 5kil, the maximal resistance 
of a single graphene sheet. The differential conductance 
dl/dV was either measured by modulating the voltage 
bias and measuring the current modulation by standard 
lock-in detection or deduced from the differential resis- 
tance obtained by applying a small ac current of typical 
amplitude 1 to lOnA superimposed to a dc current bias. 
The dc voltage drop V through the sample was then de- 
duced by integration of dV/dl = f{I)- The triangular 
shape of dl/dV = f{V) observed at high bias (above 0.15 
V, see figllj can be related to the linear dependence of the 
density of states I'iE), characteristic of the band struc- 
ture of graphene as well as of graphite at high enough 



energy [l]|: Indeed, the electronic transmission between 



the graphite sample and the underlying electrodes is low, 
so that the voltage drop occurs mainly at the contacts. 
The differential conductance can then be written as: 



dl/dV oc TLTR{av{EF+aeV) + {l-a)v{EF-{l~a)eV)) 

(1) 
where T l and T ji are the transmissions of the left and 
right contacts respectively. The parameter 1/2 < a < 1 
characterizes the asymmetry of the contact resistances 
and voltage drop, a = 1/2 corresponds to symmetrical 
contacts with V/2 voltage drop at each contact. The 
asymmetry observed between positive and negative bias 
is attributed to a combination of a slight doping of the 
sample together with some asymmetry in the transmis- 
sion of the electrodes. 

We now focus on the conductance at low voltage (be- 
low 0.12 V) which exhibits a pronounced dip at zero 
bias. This behavior is characteristic of Coulomb blockade 
through a small capacitance tunnel junction in series with 
a dissipative electromagnetic environment which can ex- 
change energy with the tunneling quasi particles on a 
scale much smaller than the charging energy. This yields 
the so called Dynamical Coulomb Blockade (DCB) [9|. 
The differential conductance data is expected to follow a 
scaling behavior as a function of bias and temperature: 



G[V) = dl/dV = T^f{eV/kBT) 



(2) 



with lim.x^of{x) = Cst and linijj^oo /(a;) = x^ , and 
the exponent z is expected to be o^R/Rq inthe case of 
Two asymmetric junctions, where R is the resistance of 
the environment and Rq ~ /i/2e^ is the resistance quan- 
tum. The data shown in figH] yields z of the order of 
0.25 ±0.05. Such a power law dependence also was found 
in the two other thicker samples, with similar exponents. 
It is thus reasonable to assume that the dissipative ohmic 
environment in the present case is constituted by the few 
top graphene layers of the graphite samples. Note also 
that a similar behavior was already observed on multiwall 
carbon nanotubes with low conductance contacts [l2l |. 

More original, as shown on fig. [3l is the bias depen- 
dence of the differential conductance measured on the 
thinest sample investigated which is 10 nm thick and 
contains thirty graphene layers. It exhibits a series of 
eight sharp dips resembling the zero bias one and nearly 
equally spaced by 20 ± 2mV. Their amplitude decreases 
with increasing voltage except for the broader dip at 50 
mV which can be decomposed into 2 overlapping nega- 
tive peaks centered around 40 and 60mV as suggested by 
the data taken at InA ac excitation. The energy scale of 
20 meV does not correspond to any simple electronic en- 
ergy scale in the sample, whose charging energy is in the 
meV range and level spacing in the lOfxeV range. On the 
other hand the lowest energy optical phonon in graphite 
(so called ZO') has an energy of 15meV [8|]. This mode, 
which emerges from the out of plane transverse acous- 
tic mode of graphene [8|, |l3| , is only present in graphite 
and corresponds to the two neighboring, non equivalent, 
planes vibrating in phase opposition along the c axis. 
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FIG. 1: Bias dependence of the differential conductance mea- 
sured on a suspended thin foil of graphite. Inset: a) transmis- 
sion electron microscopy picture of the sample, b) Side view 
taken at high resolution from which it is possible to estimate 
the number of graphene layers of the order of 30. 
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FIG. 2: Differential conductance in the vicinity of zero bias 
measured on a thick suspended graphite sample measured at 
several temperatures between 300mK (lower curve) and IK 
(upper curve) The continuous line is a power law of expo- 
nent 0.25. The data can be rescaled according to eq[2]with z 
=0.25. Full circles: temperature dependence of the zero bias 
resistance. 



This phonon mode has been observed experimentally in 
graphite by inelastic X-ray scattering 13j and scanning 
tunneling spectroscopy |14| with an energy of 15±lmey. 
The observed peak positions at multiple values of 20 
mV instead of 15 mV can be attributed to the asymme- 
try of the contacts which corresponds to the parameter 
a ~ 0.75 in eqlT] These dips are only observed on the 10 
nm thick foil and were not detected on the two thicker 



(more than 100 nm) samples. This can be understood 
considering that the conversion from electric energy (de- 
pending only on the resistance of the tunneling barriers) 
into mechanical vibrations leads to an induced vibration 
amplitude inversely proportional to the number of lay- 
ers in the graphite foil. The suspended character of the 
sample is also essential, since interaction with a substrate 
suppresses considerably the amplitude of induced vibra- 
tions as already demonstrated on carbon nanotubes [3(. 
Note that STM spectroscopy on bulk graphite [lj| also 
reveals inelastic contributions due to plasmons which are 
not detected in the present experiment. 

In order to explain the data more quantitatively we 
extended the work of Mitra et al. [7|] on the phonon 
assisted Coulomb staircase observed in the transport 
through fuUercncs molecules. The coupling between the 
ZO '-phonon mode and the electrons in the graphite sam- 
ple is described using a Holstein Hamiltonian 15| . In the 
absence of disorder this Hamiltonian Hq reads: 



Hn 



y^ e-kclck + \hw'^ c^Ckib'^ 



+ b) + hLub+b (3) 



where c^ and Ck are the fermionic electron creation and 
annihilation operators in momentum space, Ck the elec- 
tronic energy, and 6+ and b the creation and annihila- 
tion operators of the bosonic ZO'-phonon of frequency 
uj. In contrast with previous work [7[, the electronic 
energy level spacing is small compared to the phonon 
energy fko. The parameter A is the dimensionless elec- 
tron phonon coupling constant which we assume to be 
of the order of unity like in carbon nanotubes [16| . The 
coupling to the leads is then described in a dynamical 
Coulomb-blockade formalism [9[ , by a Hamiltonian of the 
form Ht = Ylik,k' Tk,k'akcl,e~'^'^ + he. Here ak are the 
electron annihilation operators in the leads, Tk^k' are the 
tunnel amplitudes and is a phase operator describing 
the electromagnetic environment of the junction. 

The Hamiltonian Hq can be diagonalized with a 

s 

In 

the limit when the charging energy of the sample is 
negligible [7| the transformed Hamiltonian reads simply 
Hg = Ylik^kCk'^k + Sw 6+6 where we have omitted the 
primes for the transformed operators. In the transformed 
basis the transfer Hamiltonian Ht is given by : 



canonical Lang-Firsov transformation: 6' = e ^6e 
H'g 



'Hae^ where S = XJ2k4ck{b - 6 



Ht = y^Tk,k'akCf.,e 

k,k' 



+ ~i4>+Hb-b+) 



he 



(4) 



This expression is obtained by expanding the product 
e~^ Htc^ under the assumption that the environment 
phase (j) commutes with the phonon operators. It shows 
that the coupling to phonons essentially changes the 
phase operator of the junction. As a result the cur- 
rent through the junction can be expressed with an ef- 
fective P{E) function describing the probability of elec- 
trons to loose an energy _E in a tunnel transition as in 
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FIG. 3: (a)Differential conductance measured at 70mK on 
the thin graphene foil depicted in fig.??, for to current exci- 
tations InA (thin line) and lOnA (bold line). Note the sharp 
dips nearly equally spaced by 20mV , except the broad dip 
centered around 50mV which can be decomposed into 2 over- 
lapping dips centered around 40 and 60mV as suggested by 
the data taken at InA ac excitation. (b) Theoretical fit for 
different phonon temperatures with the following parameters 
deduced from the geometry of the sample and the differential 
conductance data at low bias described by DCB : a = 0.75, 
Rt = 125kn, A = 200mV, e'^/C = huj/4:. The environment 
is described with a resistance R/Rq = 1. The only free ad- 
justable parameter is A = 0.7. 



usual DCB theory. Since the electromagnetic environ- 
ment and phonon operators commute, this function can 
be expressed as a convolution 



P{E) - / dE'Penv{E')Pph{E~E') 



(5) 



where Penv{E) is the probability of emitting a photon 
of energy E in the RC environment of the junction and 
Pph is the probability of emitting a phonon in the sam- 
ple. The probability distribution Pph can be obtained by 
noticing that the corresponding phase operator iA(6— 6"*") 
is analogous to that of an electromagnetic LC circuit with 
resonant frequency -^= = cj [9| (this result can also be 
obtained directly by tracing out the phonon degrees of 
freedom in eql4|): 



Pph{E) 



_ -A2coth(i 



y 6{E — khujo)e 



kphuj/2 



Ik{ 



A2 



sinh(/3fiw/2) 
(6) 



where (3 is the inverse of the phonon temperature 
Tph- Using Eqs. ([5] and [6|) and standard expres- 
sions for current as a function of P{E) [91] it is 
possible to compute the current through the sam- 
ple. For example in the case of symmetric contacts 



{a — 1/2) and constant density of states, dl/dV = 
^ J dE {P{eV/2 - E) + P{-eV/2 - E)) where Rt is 
the tunnel resistance of each contact. Note the similar- 
ity with the DCB in a tunnel junction in series with a 
LC resonator [10| . The case of an energy dependent den- 
sity of states , and of asymmetric contacts can readily be 
included by straightforward generalization of this expres- 
sion according to eq[TJ For comparison with the experi- 
ments we assume that the density of states of graphite is 
of the form: p{E) = uq{1 + '-^). This formula is exact for 
bilayer graphene with A sa AOOmeV 17| and \E\ < A. 
We also take the values of Rt, environment resistance 
R and charging energy deduced from the geometry and 
conductance data at low bias. The only free adjustable 
parameter is A = 0.7. As shown on figOthe agreement 
with experimental data is only qualitative especially at 
zero temperature , where theory predicts zero conduc- 
tance at zero bias which is not observed experimentally 
. Better agreement is found when a finite phonon tem- 
perature is introduced. Since the sample is suspended 
the populations of phonons are supposed to be strongly 
out of equilibrium. We have tentatively introduced a 
phonon temperature increasing linearly with bias which 
leads to better agreement with experimental data. How- 
ever, surprisingly, the best fit is obtained by imposing 
a bias-independent phonon temperature which does not 
seem a priori very physical. We have also in our model 
neglected the electronic temperature which is expected 
to be in the K range. Moreover we have not included the 
expected broadening of the phonon modes (even at very 
low temperature) due to the strong coupling to electrons. 
This may explain why the best fit corresponds to a finite 
phonon temperature. 

In the following we discuss the magnetic field depen- 
dence of the dips in the differential conductance ( fig. 
|4|) . They vary both in amplitude and position between 
and 5 T . Whereas the first two dips and the forth one 
are shifted toward lower frequency with increasing mag- 
netic field (as seen both at positive and negative bias) 
the third peak is shifted toward high frequency (we con- 
sider here the negative bias data since the second and 
third dips at positive bias can barely be resolved as dis- 
cussed above). The amplitude of the dips (see inset in 
figHl) may decrease or increase with magnetic field (peak 
2 and 3), or vary in a non monotonic way (peak and 
1). The relative shifts of the dips with magnetic field 
by an typical amount of 5 to 20% are of the same order 
of magnitude as the relative variations of their magni- 
tudes. We attribute these effects to the field dependent 
density of states of the graphite foil in the field range 
where Shubnikov de Haas oscillations just start to show 
up. Even if these observations are not yet understood 
in detail they indicate a strong electron phonon coupling 
and justify the value of A = 0.7 since a relationship such 
as Au;{B)/u;{B = 0) = \^lov{Ef)AG{B)/G{B = 0) 
between the typical phonon frequency magnetic depen- 
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FIG. 4: Evolution of the differential conductance dips with 
magnetic field. Insets: field dependence of the amplitude and 
position of the dips labeled 1,2 and 3. 



dence and magneto conductance, is expected to hold 
[18J . High electron phonon coupling has already been 
reported in graphene concerning in plane optical modes 
[l9|. In the present case however, strong electron phonon 
coupling between transverse ZO' vibrations and strongly 
anisotropic transport in the thin graphite layer is not 
straight-forward but can be understood if the electrical 
contacts between the two electrodes and the graphite foil 
take place through distinct graphene mono-layers which 
is highly probable. 

In conclusion we have shown evidence of differential 
conductance sharp dips in a suspended thin layer of 
graphite with 30 graphene foils. These dips can be in- 
terpreted within a simple model of dynamical Coulomb 
blockade with an environment strongly coupled to the 
lowest energy optical phonon mode ZO' of graphite. The 
magnetic field dependence of the effect corroborates this 
interpretation. 
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